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Abstract 

The conditions for the existence of heterochnic connections between the transverse 
modes of a CO2 laser whose setup has a perfect cyhndrical symmetry are discussed 
by symmetry arguments for the cases of three, four and five interacting modes. 
Exphcit conditions for the parameters are derived, which can guide observation of 
such phenomena. 
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1 Introduction 



Recent evidence of space-time complexity in lasers suggests the possibility of using 
them as "test-benches" for space time chaos theories [Arecchi et ai, 1990,1991,1999] 
ID'Angelo et a/., 19921 riredicce et a/., 1991l[nreen et a/., 1990| . Most lasers used in 



applications operate with a single Gaussian transverse mode. However, several labo- 
ratory configurations dealt with many transverse modes simultaneously active within 
the optical cavity. The nonlinear coupling between them gives rise to a complicated 
dynamics, and there is hope that a careful study of the transition from a situation 
where only few modes are active to a situation with many active modes will shed 
light on the mechanisms of transition to space time chaotic states. High accuracy 
and short times involved in laser phenomena give to lasers a relevant role in the study 
of such kind of processes. Presently, there is a wide spectrum of experimental fea- 
tures observed in transverse laser patterns waiting for theoretical explanation. These 
range from configurations in which only a few modes are active to ones in which the 
interaction among many modes leads to a loss of space and time correlation. An 
example in the first group appears in a laser system with photorefractive materials 
[Arecchi et ai, 1990| . There, the competition between different transverse modes 
can give rise to periodic or chaotic alternation between them. This phenomenon has 
been theoretically explained [Arecchi et ai, 1992| by symmetry arguments point- 
ing to the existence of heteroclinic connections involving the first three transverse 
modes, with the assumption that a subcritical bifurcation is responsible for the birth 
of the fundamental Gaussian mode. 

In this paper we consider explicitely a laser system in which the transverse modes 
are all born supercritically. Our main aim is to discuss the possibility for obtaining 
heteroclinic cycles when three, four or five modes are active in a CO2 laser. The 
paper is organized as follows. In section 2 the existence of heteroclinic connections 
in dynamical systems with symmetry is reviewed. In section 3, the particular laser 
system under study is described. The analysis of the mode equations for this system 
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for a perfect symmetry is done in section 4. The conditions that the coefficients of 
those equations must obey in order to have some kind of heterochnic connection in 
a CO2 laser system are obtained. The main conclusions and a short discussion on 
possible experiments are gathered in section 5. 



2 Symmetry Induced Heteroclinic Cycles 



Let Pi, P2,---, Pn be fixed points of a system of ODEs. A heteroclinic cycle (HC) is 
a set of trajectories Qi{t),Q2{t), Qnif) such that 

(1) 
(2) 
(3) 

Such a connection between invariant sets is non generic for systems of ODEs. But 
if the system has some symmetries, it has been shown that HCs can be structurally 
stable Field, 19801 Melbourne et ai, 1989| . The basic idea is that in the presence 
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of symmetries, there are subspaces invariant under the fiow. The role of these fiow 
invariant subspaces in building HCs is determinant. Take, for instance, a particular 
example relevant for the system discussed in this paper. Let 

x = f{x) (4) 

be a three dimensional dynamical system with x G R^, Hi, i = 1,2, 3, three invariant 
planes and Qi(t), i = 1,2,3, three trajectories connecting the invariant planes (see 
Figure 1). Suppose that the saddle fixed points Pi, i = 1,2,3, are the asymptotic 
limits of the trajectories Qi(t) such as indicated in the relations ((DEI). Thus, the un- 
stable set of the fixed point Pi is contained in the invariant plane Llj and its stable set 
is contained in the invariant plane nj_i. Therefore, we have a heteroclinic connec- 
tion between the points Pi, i = 1,2,3, forced, in this case, by the symmetries of the 
system and the boundedness of the trajectories. Let us assume that the planes are 
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induced by subgroups of the symmetry group of the problem. If the system is per- 
turbed with a term equivariant under the same group, it will have the same invariant 
planes, and, as saddle-sink connections are structurally stable in R^, the cycle will 
persist. The situation is quite different if the system is perturbed with a non equiv- 
ariant term. In that case, if the flow invariant planes are destroyed, one can break 
the heteroclinic cycle and induce the appearance of periodic orbits, giving rise to 
the phenomenon of periodic alternation, as it has been recently shown, for instance, 
in a rotating convection system |Lin fc Vivancos, 2002| . We are interested in the 
possibility of a periodic alternation among modes with different angular momenta 
for a CO2 laser in a Fabry-Perot cavity. In a perfect setup, the boundary conditions 
of this system (cylindrical tube) imposes an 0(2) symmetry. It has been shown that 
slight imperfections in the system (a non perfect coaxial discharge, slight disalign- 
ments of the mirrors, etc.) can break the rotational component of 0(2), leaving just 
the Z2 component (reflections) |D'Angelo et al, 1992| [Lopez- Ruiz et al., 1993,1994]. 
The existence of heteroclinic cycles connecting invariant sets of a 0(2) equivariant 
system was studied by Melbourne et al. [Melbourne et al., 1989| . Those cycles were 
found in the truncated normal forms of Hopf-Hopf mode interactions. In the fol- 
lowing section the normal form describing the interactions between modes in the 
OO2 laser system is presented and the relationships among the coefficients to have 
heteroclinic cycles are determined. 



3 The System 



The system under study is a OO2 laser in a Fabry-Perot cavity [Solari fc Gilmore, 1990| . 
The cavity has a perfectly reflecting plane mirror at one end (z=0) and a curved mir- 
ror with partial reflectivity at the other (z=-l). Physically, the effective curvature 
of this mirror can be modified inserting a passive optical device. As the pumping 
strength is changed, different modes will be excited. The resulting dynamics can 
be studied by reducing the Maxwell-Bloch equations to evolution equations for the 
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amplitudes of the active modes. This reduction procedure is made in two steps. 
First, the electric field, the polarization and the population inversion are expanded 
in terms of left and right longitudinal modes 0"*", 



p = z(p)e-*^[(0+rn„ + (rrn_„](52+ ^2)1/2^ 
D = (L')(L'o+E[^n(0^r+^-n(rr]), 



(5) 
(6) 
(7) 



where {E),{P),{D) stand for scale factors, x is the rate of cavity losses, s is a 
constant that depends on the curvature of the mirrors and z indicates the axial 
direction. Evolution equations for the slowly varying amplitudes in that expansion 
are derived. They read 

1 r TT+ I.-. . A nr^n r-rrn 

(8) 
(9) 
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(10) 

(11) 



where 7 is the rate of decay of the population inversion, /? is the atomic rate of 
decay of the atomic polarization, K is the pumping profile, H~^~ are 3D differential 
operators and I2 is the 2D identity matrix. The trivial solution of these equations 
is 



= 0, 

n+- = 0, 

D+-2 = 0, 



(12) 
(13) 
(14) 



and 



Do^K. 



(15) 
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For certain values of the parameters, different cavity modes are born from a Hopf 
bifurcation of this trivial solution. The electric field amplitude can be expanded as 
a combination of those cavity mode amplitudes: 



(16) 



Introducing this expansion into the Maxwell-Bloch equations (jS HTlj) one arrives to 
a set of evolution equations for the cavity mode amplitudes: 



Zo, = Laf^Z/B + Maf^ul3Zf,zlZi3 + k.O.t. 



The coefficients in these equations can be explicitly computed: 
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(18) 
(19) 
(20) 

(21) 
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s"^ + z"^ 

with Via the slow temporal frequency of the empty cavity modes. 

In the following section the relationships between the coefficients of these equations 

to have a HC are derived. 



4 Solutions of the Amplitude Equations 

Now the conditions for the existence of heteroclinic solutions in the CO2 laser system 
are discussed. If the Gaussian mode interacts with the modes of angular momenta 
±1, HCs are not possible. This result is generic for 0(2) equivariant systems when 
all modes are born supercritically from the origin. Then one must go one step further 
considering four modes (of angular momentum ±1, ±2) interacting in the cavity. 
It is interesting to analyze this case because this is the simplest situation in which 
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heteroclinic cycles can be built generically for 0(2) equivariant systems. However, 
due to the specific form of the physical equations we prove that HC are not possible 
in a CO2 laser in this case. The next step is to consider five modes with angular 
momenta 0, ±1, ±2. The restrictions among the coefficients of the normal form to 
have a HC are determined. We separate the discussion of each case in separated 
subsections. 

4.1 Three Mode Interaction 

The equations describing the interaction between the modes with angular momenta 
0, ±1 near the Hopf-Hopf bifurcation (the parameter value for which the primary 
branches in which the modes are born cross each other) are: 

ii = \^z^-{A{z^zI + 2z2zI) + DzqzI)zi, (22) 
Z2 = XiZ2-{A{2z,zl + Z2z;) + Dzoz*)z2, (23) 
zo = XqZo - iE{zizl + Z2ZI) + BzqzDzq. (24) 

The coefficients A, B, E, D have been explicitly computed in the previous section. 
Setting Zi = the amplitude equations read 

pi = X[p,-{A^{pI + 2pI) + D^pI)p,, (25) 
P2 = X[p2-{A^{2pI + pI) + D^pI)p2, (26) 
po = yoPo-{E^-{pl + pl) + B^-pl)po. (27) 

The codimension two bifurcation diagram, (Aq, A^^), for these equations can be found 
in figure 1 in Ref. |L6pez-Ruiz et al, 1993| . In that figure, projections into the 
amplitude subspaces are shown for different regions of the parameter space. Only 
the real part of the amplitudes must be considered because the dynamics of the 
phases is trivial. Notice that the dimension of the amplitude space is three, that 
there are three invariant planes (each defined by Pi=0) and one fixed point in each of 
the three invariant lines. But heteroclinic connections cannot be established between 
these fixed points due to the residual symmetry that remains in the plane [pi, P2) 
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after projecting in the amplitude space. This symmetry is given by the invariance 
of the equations ()25II27|) under the variable change: {po,pi,p2) {po,P2,Pi)- If 
the flow is projected on the plane (pi,p2) the dynamics derives towards one of the 
modes, pi or p2 (or pi = p2 if these amplitudes were the same at the beginning). 
In this situation one of theses modes interacts with po and there is no possibility to 
construct a HC because the origin is a repulsor and there are only two degrees of 
freedom (fig. 2). 

If one assumes a subcritical bifurcation for birth of the fundamental mode (Aq < 
in Eqs. ()25II27|) . this allows to obtain heteroclinic connections between the three 
transverse modes considered here [Arecchi et ai, 1992| . Since however we are inter- 
ested to study the CO2 laser, this cannot be assumed; on the contrary, all modes 
are born supercritically from the origin and thus all A's have a positive real part. 



4.2 Four Mode Interaction 



Once excluded the interaction among three modes that bifurcate supercritically we 
analyse the possibility of a heteroclinic connection among four modes with angular 
momenta ±1, ±2. 

According to Eq. (fTTjl . the normal form describing the interaction between four 
modes of angular momenta ±1, ±2 are 

ii = Xizi - iA{zizl + 2z2zl) + B{z^zl + zazI))zx, (28) 

Z2 = X1Z2 - iA{2ZiZl + Z2Z;) + B{Z3Z; + Z4ZI))Z2, (29) 

is = X2Z3-{D{z3z; + 2z^zl) + C{z,zl + Z2z;))z3, (30) 

Z4 = X1Z4 - {D{2Z3Z^ + Z4ZI) + (7(^3^3 + Z4ZI))Z4. (31) 

The behaviour of this system can be analyzed setting z^ = piC^'^^ and writing for the 
real part of the amphtudes the following equations, 

pi = X[p^-{A^{pl + 2pl) + B^{pl + pl))p^, (32) 



P2 = Xlp2-{A^{2pl + pl) + B^^ipl + pl))p2, (33) 
P3 = XIps-{C^{pI + 2pI) + D^{pI + pI))p,, (34) 

P4 = xip,-{c^{2pI + pI) + d^'{pI + pI))p,. (35) 

Let us observe the residual symmetry in the amphtude space: (pi, P2, Ps, P4) 
(P2! P15 P47 Ps)- The same hue of reasoning of the preceding subsection can be fol- 
lowed: the symmetry in the plane (pi, P2) (p2, Pi) reduces the dynamics to one of 
these modes. But the interaction of this mode with (ps, p^) still presents the residual 
symmetry (p^, ps, P4) — > (pj, p4, ps). Returning to the case of three interacting modes 
discussed in section 14.11 we also conclude that there is no way to construct a hetero- 
clinic cycle among these modes. Let us remark that symmetry in these equations is 
crucial in avoiding a heteroclinic connection. Thus the dynamics allows to connect 
the plane (pi,p2) with the plane (p3,P4), but it dies there without any possibility to 
return to the first level. Therefore, at least, another mode must be active in order 
to allow the system to come back to the initial level (in some sense a 'pump' mode). 
In this form a heteroclinic cycle could be obtained with a minimun of five modes. 

4.3 Five Mode Interaction 

Thus, we analyse the interaction between modes with angular momenta 0, ±1, ±2 
in a system with perfect 0(2) symmetry. In this case symmetries could induce 
heteroclinic cycles. The bifurcation equations for this case read as follow: 





— XqZq — 
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+ CzqzI)zq, 
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(39) 


z^ 


= X2Z3 — 


{G{2z^zl + Zizl) 


+ I{zizl + ^2-22) 
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(40) 



A complete study of this set of equations is beyond the scope of this paper. We 
focus our attention to the possibility of obtaining the situation displayed in figure 
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3. Let us first find the restrictions in the coefiicients to have a connection between 
the rotating waves (RW) of two interacting modes (Fig. 4): 



ii = XiZi - {Aizizl + A[z2Z2)zi, 

Z2 = X2Z2 — {A2Z2Z2 + A'2Zizl)z2. 

After projecting in the amphtudes space {zi — pie^'^') these lead to: 

Pi = A>i - AirpI - A'irpIpi, 
P2 = A2P2 - A2rpI - A2rpIp2. 

To obtain a connection as in Fig. 4 the unstable manifold of RWi = 



(41) 
(42) 

(43) 
(44) 



0) must 

coincide with the stable manifold of RW2 = (0, -j^)- This is reached by imposing 
that the eigenvalues of the Jacobian matrix arc one positive (in RWi) and another 
negative (in RW2)- As the fiow is bounded and there exist no fixed point in the 
interior of the plane (pi, P2), the unstable manifold of RWi is the stable manifold of 
RW2. The following conditions must hold: 
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(45) 
(46) 



These considerations are now extended to the five mode interaction in order to find 
the conditions to have a global heteroclinic connection of the type (Fig. 5): 



Po plane{p3, p^) plane{pi,p2) Po 



(47) 



These conditions are: 



P0^(P3,P4) < 



(P3,P4) ^ (Pl,P2) 



(Pl,P2)^P0 < 
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Two different situations can be distinguished. The first one corresponds to A"^ < 
0, if < and < 0. In such a case, condition P5|l is automatically verified for 
the three connections and those relationships ()48II50|) reduce simply to 



1 (KU 

Let us call this case " strong" HC among transverse modes. The physical meaning of 
this case clearly appears looking at Eqs. ()36M0|) . Indeed, selecting < 0,H'^ < 
and < 0, any growth of po produces a source term in the equations for (p3,p4); 
any growth of {p^, p^) produces a source term in the {pi, P2) equations and any 
growth of {pi, P2) gives rise to a source term in the po equation. So the cycle can 
close itself. 



Guided by such a suggestion, the values of nonlinear coefficients of Eqs. 
are calculated directly from Maxwell-Bloch equations following the procedure de- 
scribed in Sec. El and more explicit in the appendix of Ref. |L6pez-Ruiz et ai, 1994| , 
with the physical parameters of a CO2 laser. Calculations show that any choice of 
those physical parameters leading to a negative value of leads also to a negative 
value of A[j^. Therefore, the conditions of a "strong" HC in a CO2 laser cannot be 
fulfilled. 

The second possibility will be called "weak" heteroclinic connection. This corre- 
sponds to a positive real part for all nonlinear coefficients of Eqs. ()36H40j) . In such 
a case, the set of conditions ()48ll5()j) written above are equivalent to: 

M = Mino Mini Min2 > 1, (52) 

where 

Mmo = Mm(^,^), (53) 

Mmi = Mm(^,-^), (54) 

Mtn2 = Mtn{ — ,—). (55) 

Calculations with realistic parameters under usual conditions in a CO2 laser with a 
perfect 0(2) symmetry show numerical evidence that M has a maximum value equal 
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to 1. Nevertheless, symmetry arguments do not exclude the possibility of finding 
heteroclinic cycles in modified setups, provided that condition (j52|) be satisfied. 

5 Conclusions and discussion 

The conditions for the existence of heteroclinic cycles in a CO2 laser with a cylin- 
drical symmetry have been obtained. From general symmetry arguments we proved 
that this kind of cycles are not possible when there are only three or four modes 
-arising from a supercritical bifurcation- interacting into the cavity. 
At least five modes must interact to produce that kind of cycles. In this case the 
relationships that the coefficients must satisfy for having heteroclinic connections 
are found. Calculations for CO2 laser systems show that those requirements can not 
be fulfilled under usual conditions. 

Some modification in the nonlinear coupling of the Maxwell-Bloch equations is re- 
quired, such as introducing a nonlinear passive device into the cavity that modifies 
the nonlinear coefficients. 
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Figure Captions 



1. Hipotetic sketch of a possible heteroclinic connection {Qi,Q2,Q3} for a 3D- 
symmetric system of ODEs. Eacli line Qi{t) connects the fixed points (Pj, Pj+i) and 
lies on the correspondent invariant plane flj. 

2. 2D-projections of the flow in the amplitude space for the interaction in a CO2 
laser of modes and ±1. 

3. A fifth mode, po, connecting the two planes of symmetry, (pi,p2) — {P3,P4.), 
of equations ()32II35|) is needed in order to build a heteroclinic connection. 

4. A generic heteroclinic connection between the rotating waves of two interact- 
ing modes. 

5. Heteroclinic connection established with five interacting modes in a CO2 
laser: po plane {ps, p^) plane {pi, P2) po- 
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